Fluctuation symmetries for work and heat 



o 
o 

(N 
X> 

PL, 
O 



Marco Baiesi, Tim Jacobs0 Christian Maes, and Nikos S. Skantzos 
Instituut voor Theoretische Fysica, Celestijnenlaan 200D, K. U.Leuven, B-3001 Leuven, Belgium 

We consider a particle dragged through a medium at constant temperature as described by a 
Langevin equation with a time-dependent potential. The time-dependence is specified by an ex- 
ternal protocol. We give conditions on potential and protocol under which the dissipative work 
satisfies an exact symmetry in its fluctuations for all times. We also present counter examples to 
that exact fluctuation symmetry when our conditions are not satisfied. Finally, we consider the 
dissipated heat which differs from the work by a temporal boundary term. We explain why there 
is a correction to the standard fluctuation theorem due to the unboundedness of that temporal 
boundary. However, the corrected fluctuation symmetry has again a general validity. 
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I. INTRODUCTION 

Recent years have seen an explosion of results and 
discussions on a particular symmetry in the fluctua- 
tions of various dissipation functions. While started 
in the context of smooth dynamical systems and of 
thermostating algorithms and simulations p], @, soon 
the symmetry was judged relevant in the construction 
of nonequilibrium statistical mechanics. Moreover, a 
unifying scheme was developed under which the various 
fluctuation symmetries were found to be the result 
of a common feature. The basic idea is there that 
a dissipation function for a physical model can be 
identified with the source of time-symmetry breaking 
in the statistical distribution of system histories, see 
e -g- El tl f° r more details. That dissipation function is 
mostly related to the variable entropy production but, 
depending on the particular realization, can also refer to 
heat dissipation or to dissipative work. For a given ef- 
fective model, one of course always needs to check again 
that basic relation between time-reversal and dissipation. 

In the present paper, we look at a particle's posi- 
tion xt that undergoes a Langevin evolution for a 
time-dependent potential Ut- Because of that time- 
dependence, which is externally monitored, work W is 
done on the particle. At the same time, some of it flows 
as heat Q to the surrounding medium, checking the con- 
servation of energy AU — U T (x T ) — Uq(xq) = W — Q for 
the evolution during a time interval t ^ r. Both 
W and Q are fluctuating quantities and they are path- 
dependent. Our main result concerns a symmetry in the 
fluctuations of W. We give conditions on the potential 
and on its time-dependence U t under which the fluctua- 
tion symmetry for W is exact, i.e., that for all times r, 
without further approximation, 



Prob po (W dis = 
Probp (W dis = - 



exp (/3w) 



(1) 



Here, th e particle's position is initially distributed with 
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Pq ~ exp(— f3Uo) according to thermal equilibrium at 
inverse temperature f3. The notation W dls refers to the 
dissipated work (|10|l which equals the work W up to 
a difference in free energies. If the evolution would be 
reversible, then W dls — 0. In general, and confirming 
the second law, we have (W dis ) > but (JTJ also takes 
into account the trajectories where W d%s < 0. The exact 
fluctuation symmetry tells us that such "negative 
dissipative work" trajectories are exponentially damped. 

Since the heat Q differs from the dissipated work only 
by a temporal boundary term A, Q — W %s — A where 
A = A(t; xq, x t ) is non-extensive in time r, one could 
perhaps expect that Q satisfies the standard fluctuation 
symmetry, i.e., that the same as in JIJ is true after taking 
the logarithm and letting r f +00: 



1 Prob po (Q = Q T ) 



lim — log 
tT+oo r Prob 



(0 



-QT) 



= /3q 



(2) 



for q the heat per unit time. Interestingly, that is not 
what always happens, see 0. We will explain how 
the unboundedness of the potential U T can change 
the symmetry |J2J. For small enough q (basically, for 
< qr < (W dis )) the relation @ remains intact but for 
all large enough q the lefthand-side of © saturates and 
is constant. 

In what follows, we discuss the symmetry relations Q 
and @ in mathematical detail. In particular, we give 
near to optimal conditions on potential and protocol for 
which (JTJ holds. Before, that was shown only for the case 
of a harmonic potential where the minimum of the po- 
tential is moved with a fixed speed via an explicit calcu- 
lation H . There it was also found that (|2J) can be broken 
and the modification was explicitly calculated. Here we 
will argue for more general protocols and potentials to 
give estimates about the range of validity of (J2J. The 
main point is to understand when and how terms, non- 
extensive in the time r, can still contribute to the large 
deviations of the heat Q. 
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II. MODEL AND RESULTS 



In the present paper we apply the general scheme and 
algorithm of 0, Q to a model that has previously and re- 
cently been considered by a number of groups [a, IS @ • 
We find optimal conditions on potential and protocol un- 
der which the dissipated work satisfies an exact fluctua- 
tion symmetry, i.e., one that is valid for all times. The 
heat differs from that dissipated work via a temporal 
boundary term and also satisfies some general fluctuation 
symmetry asymptotically in time. Because the potential 
is unbounded, that last symmetry is not the same as in 
the standard steady state fluctuation theorem. Below we 
give more details. 



A. Model 

We consider a family of one-dimensional potentials 
Ut,t £ [0, t], as parameterized via a deterministic pro- 
tocol j t - U t {x) = C/(x,7t), with x,7t £ R. The corre- 
sponding equilibria at inverse temperature (3 are 



Pt(x) 



-pu t (x) 



(3) 



/+oo 
dx e- 0Ut{x) 
-oo 



with Helmholtz free energy F t . The time-dependence in 
7t is supposed to be given and can be quite arbitrary; of 
course the partition function Z t must be finite. The dy- 
namics is now specified by a Langevin-Ito-type equation 



dxt 



dx 



{x t ) dt + 



dbi 



(4) 



where dbt is standard white noise. Such dynamics have 
been considered before in a wide variety of contexts 
but for fluctuation theorems the emphasis has been on 
the Gaussian case. An experimental realization [j| of 
that dynamics was theoretically investigated by [j| , who 
started from (@J with 



U t (x) 



(x - vt) 2 



(5) 



A more general analysis for driven harmonic diffusive 
systems was given in Quite recently in [jj further 
experiments were considered for more general potentials 
and protocols. 

In the present paper we work with the general J2J but 
we sometimes restrict ourselves to the physically relevant 
case of 



U t {x)=U{x- lt ) 



(6) 



for a given protocol 7 = (j t , t £ [0, r]) that marks the 
shift in a potential U as time goes on. In that case F t 



does not depend on time and the associated difference in 
free energies 



AF 



(7) 



is zero. 



The model Q defines a Markov diffusion process. We 
write 



to 



(x t ,t£ [0,r]) 



for the (random) positions of the particle. If the initial 
distribution of the position xq is given via a density p, 
then Prob p (w|7) denotes the probability density of ob- 
serving a trajectory u> under the influence of the protocol 
7, with respect to the thermal noise \/2//3 dbt- Given 
a path to and a protocol 7 we also consider their time- 
reversed versions: 



e 7 , 



: 7r-i 



(8) 



B. Problem 

The observables of interest are the work and the heat. 

The work W 7 is associated to the external agent, in 
changing the potential via the protocol 7, 



W 7 (uj) 



r * ■ dUt ( \ 



(9) 



The dissipative work can then be identified with 



W* is {oj) = W-y{uj) - AF (10) 

One has to remember here that for a reversible and 
isothermal evolution the change in free energy precisely 
equals the work W 7 done on the system. Furthermore, 
in the situation © one has AF = so that W 7 = W^ ia . 

The heat Q 7 is most easily defined via the first law of 
thermodynamics. 

AU = U T {x T ) - U {x ) = W y (Lu) - Q 7 {uj) 



Q 7 (w) 



dxt 



dU t 
dx 



>t) 



(11) 

The integral (with the "o") should be understood in the 
sense of Stratonovich; it coincides better with the usual 
intuition of integrals and it does not suffer from the 
lack of time-symmetry in the Ito integral which will be 
important for us, see also [Toj |. 

In the present paper we ask 
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1 . under what conditions the work © or (|10fl satisfies 
an exact fluctuation symmetry (EFS) (J). 

2. what are the possible corrections to the standard 
fluctuation symmetry J5J for the fluctuations of the 
heat (fTTJ) . 

So far, these questions have been theoretically investi- 
gated via explicit computation for the special case of a 
linearly dragged particle in the harmonic potential O . 
in Q; question (2) has been generally addressed in [llj . 
Experimental and numerical work (in agreement with the 
results below) was done in 0, 0, IE • 



C. Results 

We start with the exact fluctuation symmetry for the 
work @. First we consider the harmonic case U(x) — 
x 2 /2 with a general protocol jt as in |JBJ. Then we give 
a general condition under which the work satisfies an 
EFS, and we give instances under which that condition is 
satisfied. Counter examples (for which the work does not 
satisfy an EFS) will show why these conditions are close 
to optimal. We end with a discussion on the relevance of 
temporal boundary terms in the large deviations of the 
heat (|TT1) . For the proofs we refer to section llVl 



1. Work 

First look at quadratic potentials, e.g. for 



U t (x) = 



(12) 



which coincides with the potential (JSJ if j t — vt. For 
that class of quadratic potentials, as in 1|12|1 . one has a 
Gaussian distribution of the work © for all protocols 
It- 

In what follows the probability density for the (dissi- 
pated) work is denoted by Prob Po (W^ s \lu) = w) as 

a function of w € JR. This (dissipated) work W 1 
pcnds on the time r, see © - <|1U|) . 



(dis) 



de- 



Proposition II. 1 (Harmonic case). 

// the distribution of the work W y is Gaussian, then for 
a general protocol 7t in © ■' 



Prob po (W,f s (ij) = w) 



Prob po (W* 5 



-w 



exp((3w) 



(13) 



for all times 



That fluctuation symmetry is easily checked to hold also 
for quadratic potentials which are more general than 
l|12p. One could argue that any Gaussian distributed 
observable can be made to satisfy a fluctuation theorem 



by rescaling the mean and the variance. However, that is 
not what happens here: no scaling at all is required for 
the work W 7 to satisfy the exact fluctuation symmetry. 

For more general potentials, we start by specifying a 
general condition: 

Assumption: We assume that there exists an involution 
s on path-space, 8=1, with sO = Os and such that 



ProVMe 7 ) = Probp,, (sw|7) 



(14) 



The involution s relates trajectories under the protocol 
7 and its time-reversed protocol O7. The next theorem 
stipulates that the existence of s implies an exact fluctua- 
tion theorem for the work. We illustrate that assumption 
below by enumerating the cases where the assumption is 
certainly verified, see also in Section ITTT1 

Theorem II. 2 (EFS Work). 

Under the assumption (|14f> above, the dissipative work 
HI U| l satisfies an EFS: for all r > and for all w, 



Prob po (W^ s (oj) = w) 
Prob po (Wf s {u) = -w) 



exp (fiii)) 



(15) 



The assumption l|14l) can be split in some two subassump- 
tions as we now state. 

Proposition II. 3. 

Suppose either (i) that the protocol is symmetric under 
time-reversal "ft — 7-r-t = @7t> or (ii) that the protocol 
is antisymmetric 7* — 70 = 7 r — 7r-t = — ©( 7 t — 7o) an d 
that the potential U obeys © and is symmetric, U(x) = 
U(—x). Then assumption (|14|) and hence the EFS <|15|) 
are verified. 

The EFS for the harmonic case U{x) = x 2 /2 with 
constant velocity 74 = vt as in © , see 0, , is treated 
by ProDOsition lII.il but is of course also a special case of 
Proposition III. 31 



We will see further in Section Till Al how our conditions 
are in fact optimal. We can however already observe here 
how some symmetry of the protocol must enter when 
dealing with an arbitrary potential. Consider indeed, if 
only formally, U(x) = x,x > with a wall U(x) = +00 
for x < in JJjJ. We can then safely assume that the tra- 
jectory satisfies Xt — jt > and © gives that the work 
Wry — W^ ls = 7 T — 7 o- Obviously this (constant) expres- 
sion never satisfies an EFS unless (and then trivially) 

7r = 70- 



2. Heat 



The heat Q 1 defined in pi|l equals the dissipative work 
W^ ls up to some temporal boundary term: 



Q 7 = W« ls + A(F - U) 



The temporal boundary A(F — U) is, modulo the factor 
/?, the change of equilibrium entropy in going from the 
equilibrium described by po to that given by p T . For the 
fluctuations of the heat we start from a situation where 
we already have the EFS (|15() for the (dissipative) work. 

We are here concerned with the situation where the po- 
tential in © is unbounded and we assume that for some 
e,v > 



U{x)>\x\ 1+e ,jt 



vt 



(16) 



at least for |x| and t sufficiently large. For the average 
work we write 



lim 



(Qj) 



+ + OC 



lim 



We further continue to assume the well-defined dynamics 
|@} with the EFS O f° r the (dissipative) work. The lat- 
ter can be summarized by introducing the rate function 
I(w) which, in logarithmic sense and asymptotically as 
t t +oo, governs 

Prob (W* is = wt) = exp[-Tl(w)} 

and assuming that I(w) > is strictly convex with 
minimum at w,I(w) = and which, from the EFS (|15|l . 
satisfies I(—w) — I(w) = f3w. Let w* be the solution of 
I'(w) — [3. In case the rate function I(w) is symmetric 
around w — w, then w* = 3w. 



Under these assumptions, we will argue in Section TlV Dl 
that the following holds true in general: 



Consider, as in @, for q ^ 
/(<?) 

Then, 



iim l lQg _Prob(Q 7 =T g ) 



t-»+oo t Prob(Q 7 = — rq) 



/(?) 



0q 

0q - I(q) 
(3w* ~I{w*) 



The antisymmetry of f(q) 
for q < 0. 



for ^ q ^ w 
for w ^ q ^ w* 
for q ^ w* 

-f(—q) fixes the behavior 



As an example, take I(w) — (3(w ~ w) 2 /Aw as is the 
case for (0, see @- Then w* = 3w and we have three 
regimes. A first linear regime where we see the usual 
symmetry (I15fl for < q ^ wJ, then a quadratic regime 
for w < q ^ Sw which saturates to a fixed value for 
q ^ 3w. Under our assumptions, we have now a general 
expression for the symmetry of the heat fluctuations, 
extending the results in |5| quite beyond the harmonic 
case (JSJ). 

A more probabilistic interpretation and a toy- calculation 
is presented in Appendix II. 




FIG. 1: Deviations from the fluctuation symmetry for the 
heat per unit time q. The function f(q) is defined in I17H . For 
small values of q, f(q) is linear so the standard fluctuation 
theorem is recovered. Between w and w* , the behavior is 
determined by the large deviation rate-function I(q) of the 
work. The function f(q) saturates for large q. 



III. EXPERIMENTS AND NUMERICAL WORK 

A. Simulations 

In the previous sections, we have given conditions on the 
protocol and on the potential for the work to follow an 
EFS (|15|) . We now argue via numerical examples that our 
sufficient conditions are also close to being necessary. To 
that aim, we have simulated the Langevin motion of the 
particle by means of an Euler-Maruyama scheme. The 
time interval [0, r] is divided into n parts dt — t /n, and 
the evolution of the system takes place via discrete states 
Xi [i = 0, 1, 2, . . . , n) connected by finite dt steps, 



Xi+l 



dU(xi,ii) 



where Bi is a random number drawn from a normal distri- 
bution, and we have set (3 = 1. The work is calculated 
through 



n-l 



W = -J2u'(xi- 7i )( 7i+1 - 7i ) 



i=0 



We consider the asymmetric potential, for a + ,a_ > 0, 

\x- lt \ a+ 

U t (x) = U(x- 7t ) 



a + 
\x-lt\ a - 



for x ^ 7 t 



(18) 



for x < "ft 



The first case examined is where the potential above is 
moved with a linear protocol "ft — t for t > 0. At t = 0, 
we generate equilibrated configurations, sampled with a 
usual Markov chain and a Metropolis criterion. First 
we chose a generic (non-harmonic) symmetric potential, 
with a + — a- — 3, for which we expect the EFS JTSJl to 
hold. That is confirmed in figure El in which we plot the 
difference 
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FIG. 2: Plot of the deviations from the EFS 1151 . for several 
values of the exponents a+ and a_ in (1181 1 . A potential which 
is dragged with constant velocity v = 1 is considered: the EFS 
is verified for the symmetric potential (here we chose a+ = 
a_ = 3), while it is not observed for asymmetric potentials. 
Parameters are r — 1 and dt = 10~ 3 . 



FIG. 3: Plot of the deviations from the EFS 1151 . for symmet- 
ric potentials [a+ = a- = a in l|18|l ] and spatially translated 
with protocol 7t = t + t 4 . We chose r = 1 and dt — 10 -4 . 
The fluctuation theorem is verified for the harmonic potential, 
while it is not valid for a symmetric potential with exponent 
a+ = a_ = 1.2. For the latter potential, a simulation with 
dt — 10~ 3 shows that numerical approximation is negligible. 



between the lefthand-side and the righthand-side in 
(|15fl . Indeed, there are no noticeable deviations from 
zero for the case of a symmetric potential. On the 
other hand, in the same figure, the results found for 
asymmetric potentials are not in agreement with the 
EFS. In that case the conditions of Proposition III. 31 are 
not verified. Note that the symmetric deviations from 
the origin found for the choices (a + = 3, a_ = 2) and 
(a+ = 2,a_ = 3) represent an indirect verification of 
the Crooks relation, see further in Section HV Al 

In figure 13 one sees again how our conditions in Proposi- 
tion HI.3l are necessary. This time we take a protocol that 
lacks the suitable temporal symmetries, like 7t — t + t 4 . 
The EFS is then not verified even for a symmetric poten- 
tial as in l|18|) with a+ = ct-{^ 2). However, as expected, 
the simulation of the special case of the harmonic poten- 
tial a+ = cy_ = 2 obeys the conclusion of Proposition 
III. II Similar conclusions are drawn from figure 0] 

B. Experiments 



-0.1 





1 

2 

y t = t 

3 

x Yt = t 


1 1 1 x - 
6 

* ^ O X 

r> x x x x 

OX X x Q _ 
X X O v 

x O o O x 

M ° X 

xe xoo 


S) ® 

n§ 9 

^xOOqOx 

U(x) = lx| 3 /3 

1,1,1,- 



12 3 

W 



FIG. 4: Plot of the deviations from the EFS (JTHJ , for the 
symmetric potential U(x) — |x| 3 /3 [a+ = a- — 3 in 1181 1 
and spatially translated with protocol jt = t 2 and jt — t 3 . 
We chose r = 1 and dt = 10 -4 . Since both protocols are 
neither symmetric or antisymmetric, the EFS is indeed not 
expected to hold. 



Previous experiments to test the fluctuation theorem 
for nonequilibrium systems included a particle dragged 
in water. In y|, Wang et al. consider a particle 
equilibrated in an optical trap and then dragged by 
the trap at constant speed relative to the surrounding 
water. The particle is micron-sized, the force is of order 
of pico-Newton and about 500 particle trajectories were 
recorded for times up to 2 seconds after initiation. The 
protocol specifies the time-dependent position of the 
trap, approximated as the position of the minimum in a 
harmonic potential with spring constant k. The external 
force exerted on the particle is thus F t (q) = —n(q — "ft)- 



The motion of 7t = vt is about rectilinear. In a second 
experiment, 0, the shape of the confining potential 
was changed. However, both are examples of harmonic 
potentials, which we have shown to be a very special 
class. 

More recently, more general situations have been inves- 
tigated, see 0,^3- E.g. a two level system was realized 
experimentally with a single defect in a diamond. When 
the system is externally driven by a laser, the dissipation 
R — (3W dis displays non-Gaussian fluctuations. It was 
noticed that integrated versions of the fluctuation the- 
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orem in their experiment are observed only for particu- 
lar protocols, in line with our general results about the 
symmetric protocols (Proposition III.3|l . In the more re- 
cent paper [6j , the distribution of the work performed on 
a particle was computed for a non-harmonic potential. 
Again, the time-symmetric protocol has been found to 
yield results consistent with the EFS. Note however that 
our results show that a symmetric protocol is not neces- 
sary; also the application of an antisymmetric (e.g. lin- 
ear) protocol combined with a symmetric potential pro- 
vides a verification of the EFS (see Proposition lII.3l and 
figure [21 ■ 



IV. PROOFS 

A. Exact identities (Crooks and Jarzynski 
relations) 

The proofs of the results listed in section III CI are dis- 
cussed here. The basic ingredient for approaching the 
fluctuations of dissipation functions via the method of 
time-reversal was already mentioned in the introduc- 
tion. In particular, for stochastic dynamics and especially 
those that we study here under equation QJ, the follow- 
ing relation is known as the Crooks fluctuation theorem, 
see ^jj; remember the notation around (jSJ): 

Lemma IV. 1. 



Prob PoH7) 
Prob PT (ew|67) 



exp (/? (Wj(lu) — AF) " 



(19) 



Proof. Using the Girsanov formula 0], the probability 
density Prob Po (w|7) on trajectories can be expressed in 
terms of the potential. Remember that the reference 
measure is associated to the U = Q case (pure Brown- 
ian trajectories), starting from po, 

Prob Po H7) 



exp 



exp 



~2 J dx t o—( Xt ) + ST 



PQJoj) 



ST 



(20) 



That leads to 

Prob po (cj\j) = p (xo) 
Prob PT (ew|67) p T {x T ) 



exp 



[/?Q 7 M] 



exp [pAU - 13 AF + /?Q 7 (w)] 



which is O since AU = -Q 1 + W 1 . 



□ 



From the Crooks relation (|19|) follows easily the so called 
Jarzynski relation |15|. In our context, this is the nor- 
malisation of the probability distribution, 



/ Prob PT (9^|e 7 ) \ /3AF 
Prob P0 M 7 ) L ' 



' po 



where (-) po is the expectation starting from p - Conclu- 
sion: 



-pAF 



i po 



(22) 



A more microscopic and physically inspired derivation of 
the Jarzynski relation follows in Appendix I. 



B. The harmonic potential with a general protocol 

For the harmonic potential i|12|) . all protocols 7 lead 
to an exact fluctuation symmetry for the work. The 
proof can easily be generalized to other quadratic forms 
of the potential where for example the protocol works 
multiplicatively (e.g. U t (x) = U{"f t x))- 

From the definition of work J^J, it is easy to see that 
the distribution of the work is Gaussian in the case of a 
harmonic potential. 



Proposition MI. ll The free energy difference Q satisfies 
AF = for all possible protocols 74. If the distribution 
of the work is Gaussian with mean w and variance a 2 , 
the expectation value in l|2^|l can be computed explicitly: 



1 = (e-WM) = 



exp 



(-2w 2 /3 + ct 4 ^) 



4 q2\ 



2a' 



where 



4 



dt 



d 2 U t . (dU t . , 



The ratio of time-forward and time-backward probabili- 
ties can then be computed by using 



Q(ST) = ST 



(21) 



Thus, necessarily, w — \a 2 fi. 

Finally, it is easy to check that a Gaussian random 
variable whose mean w and variance a 2 are related by 



W- A- 2 



a 2 P satisfies (TT3)) . 



□ 



Work EFS 



By applying property (|14fl of the involution s to the nu- 
merator and denominator of the Crooks relation l|19[) , we 



7 



find: 

exp [0{W 7 (uj) - AF)] 



Prob po (a>j7) 
Prob PT (ew|87) 
Prob PT (sa;|97) 
Prob po (9su| 7 ) 

exp [ - /3(W 7 (s6u;) + AF)] 



Hence, 



W 7 (s9u;) = -WJuj) + 2AF 



(23) 



From the first law combined with (|21[l one also concludes 
that We 7 (9w) = -W 7 (co). 



Theorem \II.2l Let us explicitly denote the dependence of 
the dynamics on the protocol 7 by writing Prob po ( W 7 = 
u>| 7 ) for the density of the work W 7 . By (|19|l 

Prob po (W/ 7 H=u/| 7 ) 

= e /3(w '- Ai " ) Prob PT (W 7 (9u;) = w'|9 7 ) 

By |jnj|: 

Prob PT (W 7 (6w) = <i/|9 7 ) 

= Prob (3o (iy 7 (9sw) = w'Y)) 

As a consequence, via (P3|) 

Prob, (M/ 7 M= W '| 7 ) 

= e ^ w '- AF *> Prob Po (W 7 (u)) = 2AF - u/| 7 ) 



Substituting «/ 
required. 



AF, we find the EFS JTSJl as 

□ 



Provosition MI. ,91 Suppose first a symmetric protocol 
97 = 7 and hence 7 T =70, 



C/ (x) = t/ r (x) 



po(x) = p T (x) 



with the distribution J3J). Choosing the identity 
operator as the involution s = 1, i.e., soj — ui, we find 
that l|14|l is obviously satisfied. 

For antisymmetric protocols 7 T -t = X — j t with X = 
7o + It we restrict ourselves to symmetric potentials of 
the form |JBJ. Cbserve then that 

U(X - x - 7t ) = U(-x + 7T _ t ) = C/(x - 7r _ t ) 



which, for i = 0, implies po(X — x) = p T {x). Choose 
therefore the involution s in 114fl as the flip suj = X — u>, 
in the sense that s(u))t = X — x t for to — (x t ). Then, 
by simple inspection from (|2U|I . again by using that the 
potential U is even, we get the equality Prob Po (sw|7) = 
Prob Px (w|9 7 ) of densities, as for Ijl4|l . □ 

D. Heat FT 

We give the arguments leading to (|17|) . Here we do not 
give a full proof. 

For very large r it is appropriate for our purpose to con- 
sider Q 7 /r = Wf s /r + A(F - U)/t as the sum of two 
independent random variables. That asymptotic inde- 
pendence can be argued for on the basis of mixing prop- 
erties of the Markov diffusion process 10}. We thus write 
formally, for arbitrary q, 



Prob(Q 7 = qr) = Prob (W^ s + A(F - U) = qr) 

dwe- T[I{w)+J{q - w)] (24) 



where Prob(W 7 As = wt) = exp[-r/(w)], Prob(A(F - 
U) = ut) = exp[— tJ(u)] in the usual sense of the theory 
of large deviations, as r ] +00. 



Hence, taking the logarithm of l|24|l and dividing by r | 
+00 takes us to 

h(q) = lim — logProb(Q 7 = qr) — — inf [I(w) + J(q — w)] 

T f w 

(25) 

and we want to compute f(q) = h(q) — h(—q). As 
I{w) is the rate function of the large deviations of the 
(dissipative) work, which we assume given and satisfying 
the EFS l|15|l . the only unknown is the rate function J. 

Clearly, always in the sense of large deviations, 



J(u) = — lim —log Prob Po 

rt + oo T 



( U{ X() )-U(x t ~vt) 



Here we assume again the independence for large r, this 
time between the variables U{x$) and U{x T — vt). Since 
Ut{x) > 0, if u > 0, then, by this independence, 



1 f +oa 

J(u) = - lim - log / dy Prob po (U(x T - vt) = (y - u)t) ~ (5% 

rT + oo T J u 



On the other hand, if u < 0, we have 



Now, the process x T — vt is stationary for large r: from 



J(u) = — lim — logProb po ([/(a; T — vt) = —ut) 

rt+oo T 
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equation 



d(xt — vt) 



-U'(x t - vt) dt - vdt + \j — dbt 



so that we can expect that for large r, x T — vt is 
distributed according to the Boltzmann statistics 



exp[— /3U(x T — vt) — (3v(x T 
have that for u < 0, J(u) - 



vt)}. As U(x) > \x 
-f3u. 



l+e 



, we 



Summarizing, in (|24H we can take J(u) = (3\u\. After 
all, it gives the probability of finding a huge energy 
difference A(F — U) ~ ut between the initial and 
the final state. It means that either Uq(xq) or U T (x T ) 
must be very large, and the energy has, in Boltzmann 
statistics, an exponential distribution. 



Finally, to obtain the results from section lTl C 2l one must 
still use that 

-I(w)+I(-w)=0w, I'(w)+l'(-w) = -0 

so that e.g. I'(—w) — —(3. It is then easily seen that 
h(q) = -I(-W) + /3(w + q) if q < -w, h(q) = -I(q) if 
— w ^ q ^ w* and h(q) = —I(w i ') — j3(q — w*) if q > w*. 
From these one computes f(q) = h(q) — h(—q). 



Appendix I: The basis of a Jarzynski relation 

Let r be the phase space on which we have a time- 
dependent dynamics |18| defined in terms of invertible 
transformations f t . One can think of a protocol 7 that 
changes in discrete steps so that a phase space point 
x G r flows in time t to ift.-yX € T with 

</?t, 7 = ft ■••/2/1, t = l,...,T 
For the reversed protocol 87, 

= fr-t+l ■ ■ ■ /t-1 fr 

We imagine a measure y, on the phase space T that is 
left invariant by ipt l7 : fi(<p^ B) — [J-(B) for BcT. Fur- 
thermore, r is equipped with an involution tt that also 
leaves /i invariant. We assume dynamical reversibility in 
the sense that for all t, 



ft TT = 7T f t 

As a consequence, n y^e-y ^ 



/r •••/r- 



t+l 



or 



<PT„*<Pt,&y 



Let us now divide the phase space in a finite partition T. 
It corresponds to a reduced description; each element in 
the partition is thought to reflect some manifest condition 
of the system. The entropy is defined a la Boltzmann as 

S(M) = In n{M), Met 



For example, in Hamiltonian systems one takes the Li- 
ouville measure as the invariant measure fi, and then we 
obtain the conventional Boltzmann definition S = In \ M\. 
We fix probability laws p and a on the elements of the 
partition and we specify the initial probability measure 
on r as 



M 



n(A n M) 
KM) 



KM) 



This probability measures A C T using p at the level 
of the partitions M of the reduced description and 
using the invariant measure p within each partition M . 
The reduced trajectories of the system are sequences 
lu = (M , Mi, , . . . ,M T ) where M l £ f, indicating 
subsequent moments when the phase space point was in 
the set Mi, i = 0, . . . , r. The path-space measure Pp i7 
gives the probability of trajectories when starting from 
r(p) and using protocol 7. 

The quantity of interest that measures the irreversibility 
in the dynamics on the level of T is (see also (11911 and 
1): 



R = In 



Pp^{M ,M 1 ,...,M T ) 
P^,e 7 (7rM T) 7rM r _i, . . . , ttM ) 



The point is that for every probability p and & on T, and 
for all M ,..., M T G f , 

R = S{M T ) - S{Mq) - In a(M r ) + In p{M ) (26) 

To show (|26|) we only have to look closer at the con- 
sequences of the dynamic reversibility. By using that 
p(B) = p,((p~^iTB), we have of course that 



/' 



n ft,*! nM T-t 



t=0 



t=0 



but moreover, by reversibility, the last expression equals 



fl ^r,7 Vt,e-y ^ M r- 



t=0 



which is all that is needed. 



f) <p-l t „M T . 



t=o 



As an immediate corollary, under the expectation Pp 



(27) 



A simple choice for the system and partition takes an 
isolated system where the reduced variables Mi refer to 
the energy of the system. We have still the freedom to 
choose p and a. Let us take p(Mo) — 1 where indeed Mo 
refers to the initial energy E. As final condition we let the 
system be randomly distributed on the energy shell E' . 
For these choices, in 'suggestive' notation, H27|) becomes 



In 



P El (E^E') 
P E >,e-y(E' 



S(E') - S(E) 
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Using that here AE — E 1 — E = W equals the work 
done, one thus recovers the microcanonical analogue of 
the Crooks relation (fH3)l . see also [lfij . 

The mathematically trivial identity l|27l) is the mother of 
all Jarzynski relations. The way in which it gets realized 
as for example an irreversible work-free energy relation 
depends on the specific context or example. We can also 
split the system from the environment. The reduced vari- 
ables (Mi) can for example be chosen to consist of the 
microscopic trajectory for the system and of the sequence 
of energies of the environment. For a thermal environ- 
ment at all times in equilibrium at inverse temperature 
(3, we thus get S(M r ) - S(M ) = (3Q where Q is the heat 
that flowed into the reservoir. On the other hand we can 
take p and a as equilibrium distributions, say of the weak 
coupling form 



p(M) 



Z 



h(E) 



where M = (x, E) combines the micro-state x of the 
system and the energy E of the environment, h(E) de- 
scribes the reservoir-distribution, U(x, 7) is the energy of 
the system with parameter 7. Similarly, 

HM) = - h(E) 

If we have that h(E ) ~ h(E T ), i.e., that the energy 
exchanges to the environment remain small compared to 
the dispersion of the energy distribution in the reservoir, 
we get from <|27|) in that context that 



-0Q-f3U(x T , lT )+l3U{x o , lo )\ 



which is a version of the Jarzynski relation l|22[l . 

Appendix II: Large deviations 

For the fluctuation symmetry we are interested in the 
large deviations of Q 7 /r from its average as r f +00. 
Such deviations can arise from two sources. First 
there are the large deviations of the work W 7 , which 
however, we know, satisfies an EFS. Secondly there is 
the possibility that AU also fluctuates to order r. This 
second effect is responsible for the deviations from the 
standard fluctuation relation @. After all, an energy 
is typically exponentially distributed and we can thus 
expect a competition with the fluctuations of the work. 

In order to clearly see the influence of the unboundedness 
of the temporal boundary, we consider here the simplest 
set-up in which deviations from the standard fluctuation 
symmetry can be calculated exactly. 

We consider a particle moving under the influence of a 
quadratic potential and a random force. For each time 



step i = 1, 2, . . . , t we take the work done on the particle 
Hi to be a random variable distributed according to a 
Gaussian of average m, and variance Vi . Let us also 
consider the analogue of the work 10 as the sum W T = 
(yi + . . . + y T ). By construction, the work per unit time 
w T = W t /t is again Gaussian with average w T = (mi + 
. . . + m r )/r and variance a 2 = (v\ + . . . + v t )/t 2 . If 
2w T — a 2 , then, automatically, the probability density 
function Prob(W T = hit) = Prob(uv = w) satisfies, for 
all r, 

Prob(u; T = w) w 



Prob(w r = — w) 



That is the (Gaussian) analogue of the exact fluctuation 
symmetry Ijl5|l for the work (that we here, by the 
previous construction, assume from the start). 

We now consider a new random variable (the analogue 
of the heat): 

Qr(w T ,yi,yr) = W T + T][(y T - a) 2 - (yi - b) 2 } 

where a, 6, r\ are real parameters, with density Prob(Q T = 
qr) . The aim of our toy model is to compute 

/(,)= limIln. Pr ° b( ^ = ^ 



r Prob(Q T = — qr) 



That can follow from f(q) = h(q) — h(—q) with h(q) the 
large deviation rate function of the heat: Prob(Q T = 
qr) ~ exp(rh(q)). The function h(q) is the Legendre 
transform of the generating function 



E(t) = lim-ln£; T (t) 

T T 



with 

E T (t) 



1 



(27r)^det = C 



(28) 

where, collectively, y = (w T ,yi,y T ) and y = (w T ,yx,y T ) 
represent their mean while C — C T corresponds to the 
covariance matrix of y. Doing the Gaussian integrals in 
(I28|l and taking the limit r — > 00 leads to 



E(t) 



\vt 2 +tw if t e 



-00 otherwise 



where = 1/277,1; = lim r a 2 = 2w — 2 lim T w T . 

We are now interested in evaluating the Legendre trans- 
form of the above, h(q) = — sup t [qt — E(t)]. The location 
of the supremum depends on whether (q—w) / v lies within 
or outside the interval [— 1+, As a result, h(q) becomes 
a quadratic function within the interval [— vt++w, vt++w] 
and a linear one outside. For the final result for f(q) one 
distinguishes between the following two cases depending 
on the value of W. 
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For vt+ < w 



/(<?) 



while for w < vt+ one has 



/(<?) 



2qU 



2 

2wU 



i- (g — w) 2 + qt± — \vt\ + wt± for g € [w — vt+, w + ft*] 



2 

2lUt* 



— w) 2 + qt± — \vi\ + wt± for q G [— w + vt+, w + vt+] 



for q £ [0, w — vt+] 
for q £ [W — uf* , wJ 
for g€ [W + t>£*, oo) 

for q £ [0, — w + vt±] 
for q £ [-W + vt+,11 
for q £ [ w + vt+, oo) 



(29) 



(30) 



The results of section III C 21 and of i-e. the Gaussian w — 1 and = 1 (i.e. 77 = 1/2). 
case where /3 = 1, are reproduced in (|30|1 by choosing 
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This appendix is discussed in a similar form in |l7j by 
one of the authors. 

In what follows, we will consider for convenience that 
Vi = v T = 1 although a more general scenario is of course 
possible. 



